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FOREST SNAIL FAUNAS FROM S. E. QUEENSLAND 
AND N.E. NEW SOUTH WALES (AUSTRALIA): 
PATTERNS OF LOCAL AND REGIONAL RICHNESS AND DIFFERENTIATION 
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ABSTRACT 


While regional snail faunas of low-latitude regions tend to be richer than similar areas at 
high latitude, there seems to be little difference between site (e.g. < 400 m°) richness, 
which often ranges between 20 and 40 species across a wide latitudinal range. Given this 
seeming paradox, it is surprising how little investigation has been made into how land snail 
species are packed into the tropical landscape across multiple scales. This paper ad- 
dresses this question by analyzing faunal lists from six regions spread across a 550 km 
extent in the temperate rain forests of eastern Australia. Considerable heterogeneity was 
observed both within similar habitat types within (ca. 40-60% faunal similarity) and be- 
tween (ca. 10-50%) different regions. The Charopidae constituted the largest fraction of 
the fauna (up to 70% of species), and were thus responsible for most of the variation in 
composition. While the number of species per region (ca. 30-50) was broadly comparable 
to Europe and North America, the rate of faunal turnover with distance was found to be 2- 
30 times greater. While geographic turnover in European and North American faunas is 
largely driven by large species, in eastern Australia turnover was largely driven by small, 
litter-dwelling species that tend to be unique to each region. The comparative richness of 
the eastern Australian fauna is thus largely related to evolutionary processes that have 
caused divergence at regional scales, rather than increased niche-packing or syntopic 
diversity at site scales. 
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America, Europe. 


INTRODUCTION 


Recent work in many parts of the world has 
revealed a great variety of regional patterns 
in land mollusc species richness and differen- 
tiation (Barker, 2005; Cameron, 2004; 
Cameron et al., 2000; Nekola, 2003, 2005; 
Pokryszko & Cameron, 2005; Seddon et al., 
2005; Stanisic, 1994; Stanisic & Ponder, 2004; 
Tattersfield, 1998.). These studies have tended 
to show that distance decay in similarity 
(Nekola & White, 1999) in regional faunas is 
low at high latitudes (Pokrysko & Cameron, 
2005; Nekola, 2005), presumably due to cli- 
mate fluctuations resulting in widespread ex- 
tinctions of forest faunas followed by 
recolonization from relatively few refugia. At 
lower latitudes, where such widespread extinc- 
tions did not occur, patterns are more com- 


plex and reflect the interaction of local envi- 
ronmental gradients and history with differ- 
ences in dispersal abilities and in situ 
speciation rates (Solem, 1984). These pro- 
cesses may allow some groups to display rapid 
distance decay rates while others demonstrate 
little. Thus, in semi-arid northwestern Austra- 
lia, large and relatively stenotypic camaenids 
generally show dramatic patterns of allopatric 
replacement, while non-camaenids are ubiq- 
uitous (Solem, 1988; Cameron, 1992). In 
mesic subtropical New Zealand, compositional 
differences tend to reflect turnover among 
small, leaf litter dwelling charopid and punctid 
snails (Barker, 2005) but without clear patterns 
of allopatry in congenerics. 

Local-scale community patterns are less well 
known, especially in tropical and subtropical 
environments where densities are often low, 
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leading to acute sampling problems (Cameron 
& Pokryszko, 2005). Contrary to Solem’s (1984) 
assertion, we now know that some tropical 
rainforest faunas are the richest known at 1 km? 
scales (Emberton, 1995; de Winter & Gitten- 
berger, 1998; Schilthuizen & Rutjes, 2001). At 
even smaller scales (i.e. 0.1-1 ha), patterns are 
even more obscured by sampling problems. As 
many tropical and subtropical forests are olig- 
otrophic and support only very low densities of 
snails, it is often impossible in such circum- 
stances to distinguish between genuine local 
heterogeneity and sampling effects. However, 
such heterogeneity could be significant both for 
speciation events and for the co-existence of 
closely related species (Cameron et al., 2003). 

The forests of eastern Australia have been 
extensively surveyed (Stanisic, 1994, 1997; 
Stanisic & Ponder, 2004). They constitute a com- 
plex mosaic of forest types, determined by sub- 
strate, temperature and rainfall, with a history 
of fluctuation in extent and connectedness. Par- 
alleling patterns seen in E. Africa (Tattersfield, 
1998; Seddon et al., 2005), individual sites may 
harbour up to 40 species, with considerable fau- 
nistic turnover occurring over relatively short 
distances. In this study, we attempt to relate lo- 
cal and regional variation in species richness 
and compositional turnover to broader regional 
patterns using both new surveys that attempt 
to overcome problems of sampling efficiency 
(Cameron & Pokryszko, 2005), and by reanaly- 
sis of some of Stanisic’s (1997) earlier results. 


STUDY AREAS AND HABITATS 


Four major study areas were analyzed with 
multiple samples, while two outlying sites were 
represented by single samples (Fig. 1). The 
Macleay Valley in northeastern New South 
Wales was the subject of a major survey by 
Stanisic (1997), covering a geographic extent 
of approximately 10,000 km? (maximum lin- 
ear distance of 70 km) with an altitudinal range 
extending from sea level to approximately 
1,100 m. As described and referenced by 
Stanisic (1997), this area contains a wide va- 
riety of forest types, including: 

(a) Sclerophyll Forests (SF) dominated by 
Eucalyptus species, have been historically 
exposed to frequent fires. While such for- 
ests can be further divided into Wet and 
Dry types, for these analyses, all have 
been lumped; both hold very poor faunas 
(Stanisic, 1997, and see below). 


(b) Temperate Rain Forests (TRF) generally 
occur at high altitudes and have higher 
yearly precipitation levels. While we have 
lumped Warm and Cool subcategories to- 
gether, most sites would fall in the Warm 
subcategory. 

(c) Dry Rain Forests (DRF) occur in isolated 
patches at lower altitudes, and typically 
contain Araucaria species in the canopy. 

(d) Coastal Rain Forests (CRF) supporting a 
high density of vine or liana taxa. 

(e) Limestone outcrops (LO) which frequently 
maintain either TRF or DRF vegetation in 
a SF or disturbed land matrix. These sites 
are usually small and isolated, and provide 
rocky shelter and high calcium levels 
needed by a number of species (Stanisic, 
1997). The vegetation is protected from fire 
by rock debris, and may be described as 
Vine Thicket, in which big trees are absent, 
but the vegetation is not fire-tolerant. 

In 2004, new samples were made in the 
Macleay Valley in TRF, DRF, and LO. Most of 
these sites represent repeat sampling of sites 
used in Stanisic (1997), and are distinguished 
from these prior samples by being referred to 
as “Kempsey”. 
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FIG. 1. Locations of study areas (italic upper case) 
in SE Queensland and NE New South Wales. 
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The three other major sampling areas were 
in southeastern Queensland. Samples near 
Kenilworth were made on the fringes of 
Conondale National Park in the valley and low 
hills‘ surrounding the Mary River, with an out- 
lier at Mapleton. Maximum linear extent was 
c. 10 km (Mapleton, ca. 30 km) with an altitu- 
dinal range from 120-250 m. All but one site 
were in DRF, the exception (site 3) being in 
Riverine Rain Forest (RRF) with large stran- 
oler figs, Ficus watkinsiana. Samples from 
Brisbane Forest, all collected within 10 km of 
each other, come from patches of Rain Forest 
near hill summits. All but one site (all above 
650 m) were classified as Wet Rain Forest 
(WRF), the exception (site 5) being in DRF at 
c. 450 m. Samples from Lamington National 
Park, in the Border Ranges, fall into two cat- 
egories: two at approximately 800 m in WRF, 
and two below 500 m in DRF. All sites are 
within 5 km of each other. Maximum linear 
extent (Kenilworth to Kempsey) was approxi- 
mately 550 km. Two additional samples were 
made between the most southerly of the 
Queensland areas (Lamington) and the 
Macleay valley/Kempsey. At Glenugie, the 
sample was made in DRF and Vine Thicket 
on basaltic talus slopes. At lluka, the sample 
came from a patch of CRF. 


METHODS AND MATERIALS 


Details of sampling methods used in the 
Macleay Valley survey are given in Stanisic 
(1997). Collections involved both searching by 
eye and the collection and subsequent search- 
ing of litter, but sampling intensity was not 
uniform between sites. From the sites 
sampled, we have chosen for further analysis 
69 in which the habitat is unambiguously in 
one of the five categories described above, 
and in which at least two species were re- 
corded. Sites from “woodland” or in other dis- 
turbed or mixed habitats are excluded. These 
lack many of the species found in less dis- 
turbed vegetation (Stanisic, 1997). The set of 
samples used here nevertheless includes all 
the native species found in the original sur- 
vey. 

In 2004, sites in the Macleay Valley were 
selected to be representative ofthe major habi- 
tats described above, excluding SF and CRF. 
In the other regions examined, sampling 


concentrated on the dominant rainforest cat- 
egories. The area searched varied between c. 
400 m? and 1 ha. Searches involved between 
two and four person-hours. About 10 |. of litter 
was collected and sieved at each site; coarse 
material was examined in situ, and the remain- 
der was bagged and examined in the labora- 
tory. All identifiable shells were extracted, 
identified and counted. Some large and easily 
identified species were counted in the field and 
left in situ. 

All identifications were carried out by J. 
Stanisic or D. Potter. Where morphospecies 
are recognised, but not yet formally described, 
they are given a reference number connected 
to examined and documented lots in the 
Queensland Museum. Nomenclature of de- 
scribed species follows Stanisic (1997). Lots 
constituting new locality records are retained 
by the Queensland Museum with other mate- 
rial deposited in the Museum of Natural His- 
tory, Wroctaw, and in the National Museum of 
Wales, Cardiff. 

Similarities between pairwise combinations 
of faunas have been analysed using the Nei 
index, in which the number of species in com- 
mon is divided by the geometric mean of rich- 
ness between both faunas. This statistic has 
a linear relationship to single step changes in 
similarity (Pokryszko & Cameron, 2005) un- 
like the familiar Jaccard index (Nekola & White 
1999), which is concave. For the Macleay 
Valley survey, we have also considered simi- 
larity as the proportion of common species to 
the least rich member of the pair, following 
exclusion of all species recorded only from 
single sites, as such unique species only add 
to dissimilarity. In all analyses, we attempt to 
take sampling error (Cameron & Pokryszko, 
2005) into account, chiefly by a conservative 
approach to the analysis of differences. 

We have used Whittaker’s Index (I, = S/a, 
where S = number of species recorded for the 
region and/or habitat within it, and @ = mean 
number of species per site within that cat- 
egory), and its variant (l „= S divided by the 
number of species in the richest site re- 
corded), as indicators of faunal differentiation 
between sites in the same region (Cameron 
& Pokryszko, 2004). The latter index elimi- 
nates the effect of site faunas that are merely 
impoverished versions of others, rather than 
containing distinctive species not found else- 
where. 
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TABLE 1. Basic data for the 69 sites used from the Macleay Valley survey (Stanisic, 1997). Single- 
tons, species found in only one site overall; Single/habitat, species found only in one site in that 
habitat. 50%+/habitat, species found in half or more of sites in that habitat. 


SF 
Number of Sites 17 
Number of Species 28 
Mean no. of species/site 3.64 
Standard Deviation 15 
Range of no.s of species/site 2-6 
W 7.69 
Imax | 4.67 
Number of Unique species 3 
Number of singletons 2 
Single/habitat 8 
50%+/habitat "o 0 

RESULTS 


Macleay Valley Survey 


Appendix 1 lists, by habitat, the species-level 
taxa found among the 69 sites considered. 
These include all species reported in Stanisic 
(1997). Details of site locations and faunas are 
available at the Queensland Museum, or in 
electronic form from the authors. Table 1 
shows the basic data for these sites by habi- 
tat. A total of 112 species, including five intro- 
duced species (* in Appendix), were found. 
None of the introduced species were abun- 
dant or widespread. = 

As shown by Stanisic (1997), SF sites have 
the poorest faunas, and LO the richest. The 
values of |, and l „ show that recorded SF 
and TRF forest faunas are exceptionally het- 
erogeneous, while those of LO and CRF are 
much less so. The high level of heterogeneity 
among sites is also shown by the high propor- 
tion (29 species, 26%) of species recorded 
from only one site, and by the low proportions 
of species that occur in at least 50% of the 
sites in any one habitat. Strikingly, no species 
occur in more than half of the SF sites, and 
only in LO and CRF sites are there more than 
10% of species in this category. Conversely, 
45% or more of species occur only once within 
SF or TRF sites, while 33% or fewer of spe- 
cies occurred only once in LO, DRF and CRF 
sites. Among the habitat faunas, only one spe- 
cies is unique to CRF, three to SF, while LO 
harbours 18. Only 10 species were observed 
within all five habitats. Of the 41 species miss- 
ing from LO sites, most are so infrequent (20 


TRF LO DRF CRF Total 
aa 9 15 6 69 
55 71 62 27 112 

7.95 27.9 12.73 9.5 

4.11 6.45 5.14 4.18 

MT AO D å am 

6.92 2.54 4.87 2.84 

3.23 1.78 2.82 1.93 
16 18 11 1 49 
10 9 7 1 29 
25 21 16 9 
3 21 3 9 


have only single reported occurrences) that 
their absence from limestone sites may be 
attributed to sampling error, given the small 
number (9) of observed sites. 

Table 2 shows the matrix of similarities be- 
tween habitats. The Nei index demonstrates 
that the strongest connection is between LO 
and DRF, whereas the maximum similarity in- 
dex identifies that CRF sites are most related 
to LO and DRF faunas. Both indices strongly 
indicate that SF faunas are differentiated from 
CRF and LO. 

For the overall fauna, the Charopidae are 
much more numerous (55 species, 49%) than 


TABLE 2. Similarity matrices for between-habi- 
tat comparisons in the Macleay Valley using (A) 
the Nei index, and (B) taking the number of spe- 
cies in common (excluding those recorded once 
only) as a proportion of those in the poorest habi- 
tat (ex unique). Values expressed as percent- 
ages. 


A: Nei TRF LO DRF CRF 
SF 48 38 48 47 
TRF 50 53 44 
LO 65 58 
DRF 56 
B: Ex unique TRF LO DRF CRF 
Sr 73 65 m 50 
TRF 69 69 65 
LO 78 88 
DRF 88 
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TABLE 3. Numbers and proportions of species in Charopidae in Macleay Valley sites, by habitat, by 
position in the quartiles of overall frequency, and by habitat range (number of the five habitats at left 
in which the species were recorded). Quartiles differ slightly in size because of ties in frequency. 
Note that the bottom quartile (4) consists of species recorded in only one site overall. 


Charopidae No. % Quartile 
SF 8 29 1 
TRF Zi 49 2 
LO 28 39 3 
DRF 24 39 4 
Total 55 49 Habitats 

1 
2 
3-5 


any other family (Table 3). SF has the small- 
est proportion of charopids in its fauna (8, 29%) 
while TRF has the largest (27, 49%). Nearly 
three quarters of all single-site species are 
charopids, as are two thirds of all single-habi- 
tat specialists. This single family thus accounts 
for most of the observed differences between 
habitats and among sites within habitats. 


2004 Surveys 


Appendix 2 lists the 149 species level taxa 
found by survey area. Six of the species re- 
corded are introduced (* in Appendix). Of 
these, only Lamellaxis clavulinus was recorded 
in large numbers, and then only in one site. 
Full details of sites, and of the species com- 
position of samples from each are deposited 
at the Queensland Museum and are available 


Total species Charopidae % 
30 10 33 
29 12 41 
24 8 54 
29 21 72 
Total species Charopidae % 
49 33 67 
25 10 40 
38 iz 82 


electronically. Table 4 shows the basic data 
for each survey area. Results for Kempsey are 
split by habitat, as the range of habitats 
sampled was much greater than elsewhere, 
and the LO sites are very different from all oth- 
ers. Of the categories with more than one 
sample, only in Kempsey LO do all sites meet 
the minimum requirements of at least 10 times 
as many individuals as species and at least 
200 individuals, and in Lamington and 
Brisbane virtually no site meets these criteria 
(Cameron & Pokryszko, 2005). Logarithmic 
rank/abundance curves tend to be concave in 
Brisbane and in some Kempsey habitats, with 
a tail of very rare species. They are more con- 
vex in Lamington and Kenilworth, where the 
pattern is very similar to that recorded for non- 
limestone Carpathian forests in central Europe 
(Fig. 2). Within the Kenilworth set, there is also 


TABLE 4. Basic data for each area in the 2004 survey. Ind/site, mean number of individual shells 
collected or recorded per site. For habitat abbreviations, see text. 


Kenilworth Brisbane Lamington Kempsey Kemp Kemp Glenugie Iluka 

DRF/RRF WRF/DRF WRF/DRF LO DRF TRE DRF CRF 
Sites 10 6 4 3 3 + 1 1 
Species 49 34 41 5l 36 11 18 11 
Mean/site 24.1 13.5 22.3 28.0 17.0 11 18 11 
SD 4.4 32 7.9 a7 6.2 Na WNia N/a 
Range 17-32 9-17 13-32 27-30 10-22 N/a N/a N/a 
lw 2.03 2.52 1.84 1:82 2.12 © Na N/a N/a 
los 1.38 2.00 1.48 1.70 1.64 N/a N/a N/a 
Ind/site 295 188 122 3053 136 35 433 248 
SD 162 46 51 2911 66 N/a N/a N/a 
Range 109-558 92-218 72-185 369-6141 60-185 N/a N/a N/a 
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FIG. 2. Logarithmic plot of relative species abundances (%) for the snail faunas in samples from 
Kenilworth, Queensland, and from Carpathian beech forests at Podlipowiec, S. Poland (Pokryszko 


and Cameron, unpublished). 


no significant correlation between numbers of 
species and of individuals (R = 0.22, P > 0.2, 
n = 10). Thus, with that exception, we may 
suspect that recorded site faunas are signifi- 
cantly incomplete. Nevertheless, the range of 
site richness, and the values of |, and |, are 
generally much lower than those seen within 
habitats in the Macleay survey. 

These indices are larger in Brisbane than 
elsewhere. In this area, five of the six sites 
are characterised as WRF, all above 650 m 
asl. The other site, DRF at c. 450 m, harbours 
six unique species, and is one of the two rich- 
est sites across the entire survey. Without this 
site the WRF sites of Brisbane have a mean 
of 12.8 species/site and |, and |. values of 
2.19 and 1.64 respectively, more in line with 
values for other areas. A similar pattern is evi- 
dent at Lamington, where the two lower el- 
evation DRF sites (with Araucaria) were richer 
in species and individuals than the two higher 
elevation WRF sites. Because of limited num- 
ber of samples, however, not enough statisti- 


cal power exists to test hypotheses regarding 
this pattern. At Kenilworth, one site (Kenil- 
worth 3) comes from Riverine Rain Forest, but 
shows no peculiarity in richness or native spe- 
cies composition from the remaining DRF 
sites, though it contained large numbers of 
the introduced Lamellaxis clavulinus 
(Subulinidae). 

Charopidae are the most numerous family 
in terms of species. As in the Macleay survey, 
they generally contribute to differences (rather 
than similarity) between areas, except in the 
single-site cases of Glenugie and Iluka. Omit- 
ting these two, the overwhelming majority of 
species are found in only one region, and only 
three (2.1%) are found in all four. 

Table 5 shows Nei similarity index values 
within and between regional habitat and total 
faunas. Fig. 3 shows the dendrogram of simi- 
larities between sites generated by UPGMA 
clustering. The Kenilworth faunas are the most 
uniform, with an average inter-sample simi- 
larity of 64% in the region. Strong distance 
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TABLE 5. Values of the Nei index of similarity (%) in the 2004 survey, A: on mean site-by-site com- 
parisons within and among areas and habitats; B: among the total faunas of each area. 


A: Nei, mean 


Kenilworth 
Brisbane 
Lamington 
Kempsey DRF 
Kempsey LO 
Kempsey WTRF 
Glenugie 

lluka 


B: Nei, total 


Brisbane Lamington Kempsey Glenugie 


Kenilworth 
Brisbane 
Lamington 
Kempsey 
Glenugie 


100 


NEI INDEX 
Wn 
= 


Kenilworth Brisbane Lamington 


64 32 
55 


52 34 
5 
Kenilworth 


Kempsey Kempsey Kempsey 


DRF 
20 5 
44 12 
JO 11 
40 


18 14 
17 27 
13 12 
15 

Brisbane 


Lo 


Lamington 


117046 9041416950395S40421941939472 


Iluka 


TRF  Glenugie Iluka 
6 18 8 
21 14 6 
16 6 5 
30 9 4 
20 14 17 
n/a 14 9 
n/a 22 

n/a 

A Kempsey 


W2B3415 


FIG. 3. Dendrogram of Nei index similarities for 2004 study sites (UPGMA). In Kempsey, sites 1, 4 
and 5 are from limestone outcrops; W (WayWay Forest) and B (Bob's Road) are based on very small 


samples (see text). 
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TABLE 6. Acomparison of the numbers of species found in the original Macleay Valley survey and in 
2004 for six sites (see text). Kempsey 4, 1 and 5 are limestone sites, 6 and 7 are DRF, 2 is TRF. 


Natural Arch Yessabah MtPleasant Blowhole Way Way Bob’s Rd 
Species Kempsey 4 Kempsey1 Kempsey5 Kempsey7 Kempsey6 Kempsey 2 
Macleay (JS) 33 40 32 14 22 8 
2004 30 27 ae 19 10 m 
Present both 25 25 21 12 9 6 
Aggregate 39 42 36 21 23 13 
Nei index (%) 78 76 74 74 61 64 


decay is apparent between regional faunas, 
with only 6% similarity being noted over the 
950 km extent from Kenilworth to Kempsey. 
All Queensland areas have low affinities to 
New South Wales areas (4-16%), while simi- 
larity between New South Wales sites also are 
surprisingly low (4-22%). While the dendro- 
gram also generally demonstrates this same 
pattern, it also shows that the Brisbane DRF 
habitat resembles Kenilworth DRF sites. Thus, 
all six species unique in the Brisbane area to 
this single DRF site are held in common with 
one or more Kenilworth DRF sites. 

The comparison of total study area faunas 
may provide more robust results by summing 
across multiple samples, increasing the num- 
ber of encountered microenvironments and in- 
dividuals (Nekola & White 1999). Similarity 
patterns for total faunas at this scale demon- 
strate essentially an identical pattern, with simi- 
larity falling to 14% between Kenilworth and 
Glenugie. Similarity of Kenilworth and Kempsey 
non-limestone sites fell to 6.3%. While 
Queensland sites maintained somewhat higher 
similarities (34-52%), inter-area similarity south 
from Lamington generally ranged from 12 to 
15%. The least similarity was noted between 
the faunas of Lamington and Iluka (5%). Thus, 
the New South Wales sites are equally dissimi- 
lar to each other as they are to the Queensland 
areas. This can be explored in more detail (Table 
4). With the exception of Brisbane, roughly mid- 
way between Lamington and Kenilworth, more 
than 40% of recorded species were unique to 
each region. The highest proportion of unique 
species (78%) was observed at Kempsey. While 
this may represent the presence of base-rich 
limestone habitats, geographic factors also ap- 
pear important as the frequency of unique spe- 
cies within the Kempsey non-limestone fauna 
remains significantly higher (68%) than that 
observed from all other sites. 


Macleay- Kempsey Comparisons 


Six of the seven Kempsey sites sampled in 
2004 had been visited in the earlier survey. A 
comparison of results is thus possible (Table 
6). Although areas searched do not exactly 
correspond, the differences between occa- 
sions show that even when sample sizes are 
large (Kempsey 1 produced more than 6,000 
shells), not all the fauna is necessarily de- 
tected. Nevertheless, indices of similarity for 
repeat samples at the same site are higher 
than those generally recorded between sites 
in the same habitat. It is clear that one sample, 
Kempsey 6, Way Way Forest, was seriously 
deficient in 2004, with only 60 individuals and 
10 species retrieved. This is reflected in its 
distant relationship to all other sites in Figure 
2. Using the combined data for each Kempsey/ 
Macleay site increases levels of between site 
similarities within each habitat (data not 
shown). The combined total numbers of spe- 
cies in each of these sites match the maxima 
for those habitats recorded in the earlier sur- 
vey quite closely. Overall, the 2004 survey 
recorded only 67 species in Kempsey, or 60% 
of those recorded in the larger Macleay Valley 
survey. As the two habitats not examined in 
2004, SF and CFR, hold few unique species, 
thus lower number mainly reflects the smaller 
number of samples made and the smaller pro- 
portion of the original survey area sampled. 
While the earlier Macleay survey clearly gives 
a more comprehensive view of the fauna of 
each habitat, it does not identify much differ- 
ent rates of turnover with other faunas. For 
example, similarity between DRF in the 1997 
Macleay survey and Kenilworth provided a Nei 
Index of 16.7% versus 18% for the 2004 
Kempsey data. Deficiencies in sampling are 
thus not seriously distorting similarity esti- 
mates between faunas. 
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DISCUSSION 
Reliability and Sampling Error 


It is clear that in both the 1997 Macleay Val- 
ley and 2004 Kempsey surveys, site faunas 
were incompletely inventoried (Cameron & 
Pokryszko, 2005), largely due to the low den- 
sities of individuals on non-limestone sites. As 
shown by the very high values of | |, and 
the large proportion of singleton records, this 
problem was most acute in the 1997 survey, 
with some of this error reflecting variation in 
sampling intensity between sites. In the 2004 
standardised surveys, this problem remained, 
but was less severe. Only Kenilworth and 
Kempsey LO sites were adequately sampled 
based on total observed richness and abun- 
dance (Cameron & Pokryszko, 2005). How- 
ever, comparisons with the 1997 survey in 
Macleay/Kempsey suggest that even samples 
of approximately > 6,000 shells may not com- 
pletely inventory all taxa, or that very small 
shifts in location may add more species to an 
inventory, even within the same habitat. The 
logistics of quadrupling sampling effort at each 
site are likely impossible within reasonable 
time frames. In spite of these limitations, 2004 
site data compare favourably with most other 
tropical studies (Emberton, 1995; de Winter & 
Gittenberger, 1999; Schilthuizen & Rutjes, 
2001), and are similar to those recorded by 
Tattersfield (1996) in Kakamega Forest, 
Kenya, where supporting evidence suggests 
that the aggregate fauna has been adequately 
inventoried (Cameron & Pokryszko, 2005; 
Tattersfield et al., 2006). The possible effects 
of such sampling error on diversity estimates 
are considered below where relevant. Com- 
pared to the traditional Jaccard index, the Nei 
index of similarity gives some compensation 
for differing richness in compared sites. Esti- 
mates (data not shown) based on the propor- 
tion of species held in common by the poorer 
of two sites shows the same pattern as that 
revealed by the Nei index. Since we are more 
concerned with faunal differences, than with 
richness per se, we have refrained from using 
rarefaction or species accumulation curves, as 
these give no indication of which taxa might 
have been overlooked. 


Habitat Effects and the Range of Local 
Richness 


Stanisic (1997) demonstrated in the Macleay 
Valley that SF was the least rich at site level, 


while LO was the richest. Our 2004 surveys 
replicate that finding, and allow estimation of 
maximum richness levels for each habitat stud- 
ied. At the site level, LO site maxima exceed 
40 species. That the total fauna of LO sites is 
considerably higher in the original 1997 (71 
species) and in the more limited 2004 survey 
(51 species) reflects both environmental dif- 
ferences between LO habitats and the pres- 
ence of single outcrop endemics (Stanisic, 
1997). Using data from both 1997 and 2004 
where possible, DRF off limestone have 
maxima of 22-23 taxa/site, TRF and CRF have 
maxima from 15-17, while SF has a maximum 
of 6. SF and CRF, both base-poor habitats, 
appear to simply harbour reduced versions of 
the faunas found in base-rich habitats, as they 
Support very few unique species, most of 
which are known from only single records. 
While more intensive work might raise these 
maximum richness figures slightly, it seems 
likely that greater between-site heterogeneity 
in base-poor habitats reflects a genuine patchi- 
ness in species distributions at scales larger 
than sample sites. We cannot tell at present 
whether this heterogeneity reflects small-scale 
geographical differentiation within the valley. 

In the Queensland series, evidence suggests 
that WRF supports less rich faunas as com- 
pared to DRF, especially in Brisbane and 
Lamington. For WRF, richness appears to 
range up to 25 per site, and to 35 over larger 
areas of a few km?. However, DRF can sup- 
port more than 30 species in single sites both 
in Lamington and Kenilworth, while they ap- 
pear to range up to 40-50 over larger areas. 

These local richness levels are similar to 
those found in forest faunas from tropical East 
Africa (Tattersfield et al., 2006; Seddon et al., 
2005), but the maxima are rather less than 
those found on the North Island, New Zealand 
(up to 64 in a site, perhaps 75+ in a few ha, 
Barker, 2005). While rainforests in Madagas- 
car, Cameroon and Borneo are richer at larger 
scales (60-110 species per 1 km*;, Emberton, 
1995; de Winter & Gittenberger, 1998; 
Schilthuizen € Rutjes, 2001), at the scale of 
400 m° plots accurate comparisons are not 
possible due to sampling errors (Cameron, 
2004; Cameron 8 Pokryszko, 2005). Surpris- 
ingly, the range in site richness values is re- 
markably similar (2-45 species/site) to that 
recorded in European (Pokryszko & Cameron, 
2005) and eastern North American (Nekola & 
Smith, 1999; Nekola, 2005) temperate forests. 
Thus, maximum site richness of forest faunas 
shows little global latitudinal variation, despite 


454 : STANISIC ET AL. 


great and deep taxonomic differences. This 
observation also holds more generally within 
eastern Australia, where sites holding c. 35- 
40 species can be found at widely separated 
latitudes (Stanisic, 1994; Stanisic & Ponder, 
2004). 

The pattern in all four of our study areas 
seems to conform to Waldén's (1981) conclu- 
sion that within a given forest area poor sites 
have merely reduced versions of richer ones. 
In areas where a greater altitudinal range is 
sampled, this not the case, with some spe- 
cies being restricted to higher or lower altitudes 
(Stanisic 1982, Tattersfield et al., 2006). 
Hence, the inclusion of poorer forest types in 
some but not all of our regional comparisons 
does not invalidate the estimation of faunal 
similarities and turnover. 


Faunal Turnover and Taxonomic Composition 


There is a considerable amount of faunal 
turnover over the approximately 550 km dis- 
tance between our most distant sites. Table 7 
demonstrates this by comparison with non- 
limestone regions along the USA-Canada bor- 
der (Nekola, 2005, and unpublished), the 
Sudetes-Carpathian chain in southern Poland 
(Pokryszko & Cameron, 2005, and unpub- 


lished), and the southern Rocky Mountain cor- 
dillera in the southwestern USA (Nekola, un- 
published). Except for Brisbane, rates of 
uniqueness in the local eastern Australian fau- 
nas were greater than 40%. Outside of east- 
ern Australia, only one other region, the 
Chiricahua Mountains (Arizona), exceeded this 
level. All remaining sites demonstrated unique- 
ness levels < 20%. Given the continuous na- 
ture of compositional turnover (Whittaker, 1975; 
Nekola & White, 1999), the arbitrary transect 
endpoints should harbour the most unique spe- 
cies, as was observed in eastern Australia, 
Poland, and the southern Rockies. The higher 
rates of unique species in the centre of the 
northeastern USA transect is due to the incur- 
sion of southern species into buffered regional 
climates adjacent to, or immediately downwind 
of the Great Lakes. 

The Nei endpoint similarity was also much 
lower for the eastern Australian transect (16.3) 
than for either of the mountain transects 
(southern Rockies = 41.7, Carpathians = 46.7) 
or for the northeastern USA (81.3). Assuming 
that similarity falls in a constant exponential 
fashion with interregional distance (Nekola & 
White, 1999), these values indicate that only 
190 km is required for faunistic similarity to 
fall by 50% in eastern Australia, while it will 


TABLE 7. Comparative data on turnover and faunal differences for four sets of sites on non-calcare- 
ous soils in E. Australia, NE USA, Poland, and SW USA (see text). Kempsey data exclude sites on 
limestone. The NE USA and Polish sequences run from west to east, the Australian and SE USA from 
north to south. The Nei ends column gives the value of the Nei index between the two faunas at the 
ends of each sequence, which are about 500 km distant except in the NE USA where they are 1,800 
km apart. 


E. Australia Kenilworth Brisbane Lamington Kempsey Nei, ends 
No. of species 49 34 41 38 1078 
Unique 20 4 18 26 
% Unique 40.8 TES 43.9 68.4 
NE USA Duluth, MN Copper Harbor, MI Northern Vermont Northern Maine Nei, ends 
No. of species 29 J9 29 19 Sus 
Unique 0 6 2 0 
% Unique 0 1832 6.9 0 
Poland Mysliborski Mlynowiec Podlipowiec Ustrzyki Nei, ends 
No. of species or 37 48 46 46.7 
Unique 74 4 5 9 
% Unique 18:9 10.8 10.4 19.7 
SW USA Chiricahua Mts. Sierra Blanca Mts. | Southern Sangre San Juan Mts. Nei, ends 

de Christo Mts. 
No. of species 25 18 20 23 41.7 
Unique 12 2 3 4 


% Unique 48.0 11.1 18570 17.4 
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take 400 km in the southern Rockies, 450 km 
in the Polish Carpathians, and 6,000 km in the 
northeastern USA. The eastern Australian 
fauna is thus experiencing geographic turn- 
over at 2.1-31.6 times the rate experienced 
by these other faunas, helping explain how the 
outlying Glenugie and Iluka sites demonstrate 
such high degrees of uniqueness: in eastern 
Australia the similarity of regions separated by 
only 100 km would be expected to fall by 30%. 

The Australian situation is similar, though not 
as extreme, as that reported by Tattersfield 
(1998) and Tattersfield et al. (1998) for Indian 
Ocean coastal forests and Eastern Arc Moun- 
tains in East Africa. In eastern Australia, single- 
site endemics are unusual except for isolated 
limestone outcrops (Stanisic & Ponder, 2004). 
Most species for which we have reliable infor- 
mation have linear ranges between 50 and 200 
km (Stanisic, 1990), and these are not always 
coherent. In this respect, these faunas re- 
semble those of North Island, New Zealand 
(Barker, 2005), and are very different from 
those in the semi-arid northwest of Australia, 
where large camaenid snails often have very 
restricted ranges (of the order of 1-10 km), 
and show strict patterns of within-genus allo- 
patric replacement (Solem, 1988; Cameron et 
al., 2005). Although camaenids and other 
large-shelled families (Rhytididae and 
Helicarionidae) do contribute to faunal differ- 
entiation in E. Australian forests, it is clear that 
the much smaller, litter-dwelling Charopidae 
contribute most, as they do in New Zealand. 
In Europe and North America small litter dwell- 
ing species tend to unify the composition be- 
tween regions as they have large geographical 
ranges (Pokryszko & Cameron, 2005). How- 
ever, in these regions it is generally the large 
helicoids and clausiliids (Europe) and large 
polgyrids, endodontids, and helminthoglyptids 
(North America) that differentiate faunas. 

In terms of taxonomic composition, we note 
that introduced species are only a small ele- 
ment of these faunas, and the few species in- 
volved are not usually abundant. While many 
of the faunas we have examined may be af- 
fected by human activities, they have so far 
retained their integrity as indigenous assem- 
blages. 


Pattern and Process 


Although interpretation is complicated by 
sampling errors in many studies (including this 
one) certain patterns can be discerned and 
reasonably attributed to known environmen- 


tal fluctuations that have influenced the de- 
gree of habitat connectedness. In both N.W. 
Australia and eastern Africa, climatic fluctua- 
tions within the Pliocene/Pleistocene have 
periodically isolated favourable forest habitats, 
within which stenotypic species have been 
confined and differentiated (Cameron, 1992; 
Tattersfield, 1998). Similar proliferation of al- 
lopatric sister species occurs elsewhere, for 
example in the Aegean region of Europe 
(Douris et al., 1998; Parmakelis et al., 2005). 
In eastern Australia, periodic isolation of mesic 
forest habitats during the Pleistocene 
(Kershaw, 1981, 1994; Lloyd & Kershaw, 1997; 
Harle, 1997) has constantly altered the extent 
and connectedness of forest types. While al- 
lopatric speciation can be pronounced on iso- 
lated limestone outcrops (Stanisic & Ponder, 
2004), as in New Zealand (Barker, 2005), few 
examples of strict allopatry or parapatry be- 
tween congeners are apparent, especially 
among the speciose charopids. The presence 
of co-existing congeners in eastern Australia 
thus suggests that speciation occurred within 
small regions, as has been suggested in east- 
ern Africa (Tattersfield, 1998). The presence 
of altitudinal zonation in the eastern Austra- 
lian fauna (i.e., Bellenden Ker, Stanisic, 1982) 
also suggests that some speciation has been 
driven, as in Africa (Seddon et al., 2005), by 
topographic complexity. 

When considering the very high local diver- 
sity of some New Zealand sites, also domi- 
nated by small charopids and punctids, Solem 
(1983, 1984) suggested that this diversity had 
built up slowly in a very stable environment. 
By implication, this process involved isolation 
by distance (which can be modest, given their 
low mobility), differentiation, and back 
colonisation into ancestral territory that was 
still occupied. This process implies that com- 
petitive exclusion is relatively unimportant, and 
that even nearby populations may be suffi- 
ciently isolated to allow for differentiation. With 
the evidence of competitive exclusion in land 
molluscs being meagre, it is worth noting that 
a metapopulation structure involving small, 
temporary and shifting populations allows lo- 
cal co-existence of ecological equivalents 
(Hubbell, 2001). These dynamics may be im- 
portant for some tropical and subtropical fau- 
nas that consist of low-density populations that 
are patchily distributed around temporarily 
favourable microhabitats, such as coarse 
woody debris (Cameron et al., 2003). As with 
estimates of local richness, further investiga- 
tion of such mechanisms will encounter formi- 
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dable sampling problems. While this current 
study certainly does not overcome these prob- 
lems, the data it generated do not contradict 
the expected patterns generated by these pro- 
cesses. 
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APPENDIX 1 


Species recorded by Stanisic (1997) in the Macleay Valley, ordered by overall frequency of 
occurrence, and showing the number of sites occupied in each habitat. For abbreviations, see 


text. Asterisked species are introduced. 


Species Sites 17 


Gyrocochlea planorbis 
Saladelos macquariensis 
Austrochloritis nambucca 
Rhytididae MV3 
Rhytididae MV7 

Thersites novaehollandiae 
Coenocharopa multiradiata 
Pleuropoma jana 
Helicarionidae MV5 

Nitor medioximus 
Punctidae MV? 

Rhytididae MV1 
Hedleyella falconeri 

lotula microcosmos 
Charopidae MV35 
Parmavitrina planilabris 
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TRF LO DRF CRF TOTAL 
22 g 15 6 69 
13 8 6 5 39 
5 sh 10 3 30 
IM 3 3 6 29 
9 7 4 2 26 
10 3 6 3 23 
T 4 T 1 21 
4 P 9 1 21 

1 9 7 3 20 
1 8 3 3 20 
3 3 10 3 19 
1 8 9) 1 15 
2 5 5 2 15 
13 1 1 0 15 
2 7 4 1 14 
3 0 a 2 14 
8 1 4 0 14 
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(continued) 
SF TRF LO DRF CRF TOTAL 

Species I7 22 9 15 6 69 
Helicarionidae MV1 14 
Allocharopa belli 13 
Tornatellinops jacksonensis 13 
Charopidae MV33 13 
Pupisoma circumlitum 12 
Pupillid MV1 12 
Hedleyoconcha delta 12 
*Discocharopa aperta let 
Charopidae MV30 10 
Gastrocopta bifurcata 10 
Charopidae MV22 10 
Camaenidae MV2 10 
Tornatellinops pressus 9 
Cralopa stroudensis 9 


Charopidae MV6 
Georissa laseroni 
Charopidae MV4 
Charopidae MV23 
Helicarionidae MV4 
Gastrocopta queenslandica 
Rhophodon kempseyensis 
Camaenidae MV9 
Glyptopupoides egregia 
Pygmipanda atomata 
Coneuplecta calculosa 
Camaenidae MV11 
Elasmias wakefieldiae 
Triboniophora graeffei 
Rhytididae MV2 
Charopidae MV25 
Charopidae MV38 
Helicarionidae MV2 
Camaenidae MV4 
Camaenidae MV7 
Cystopeltidae MV1 
Charopidae MV7 
Charopidae MV24 
Charopidae MV27 
Charopidae MV31 
Charopidae MV46 
Posorites conscendens 
Camaenidae MV8 
Camaenidae MV12 
Charopidae MV18 
Charopidae MV2 
Charopidae MV16 
Charopidae MV36 
Helicarionidae MV3 
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(continued) 


Species 


Camaenidae MV5 
Camaenidae MV6 
Camaenidae MV10 
Gastrocopta pediculus 
*Paralaoma caputspinulae 
Rhytididae MV4 
Ngairea dorrigoensis 
Letomola contortus 
Charopidae MV1 
Charopidae MV9 
Charopidae MV13 
Charopidae MV26 
Charopidae MV34 
Charopidae MV40 
Charopidae MV44 
Liardetia scandens 
Austrochloritis brevipila 
Camaenidae MV3 
*Bradybaena similaris 
Cylindrovertilla kingi 
Pupillid MV2 
*Lamellaxis clavulinus 
Succineidae MV 1 
Rhytididae MV6 


Coenocharopa yessabahensis 


Egilomen lirata 
Charopidae MV3 
Charopidae MV5 
Charopidae MV8 
Charopidae MV10 
Charopidae MV11 
Charopidae MV12 
Charopidae MV14 
Charopidae MV15 
Charopidae MV19 
Charopidae MV20 
Charopidae MV28 
Charopidae MV29 
Charopidae MV32 
Charopidae MV37 
Charopidae MV39 
Charopidae MV41 
Charopidae MV42 
Charopidae MV43 
Charopidae MV45 
Papuexul bidwilli 
Camaenidae MV1 
*Zonitoides arboreus 


Total Species 


Sites 
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Occurrence of species in each area sampled in 2004. Asterisked species are introduced. 


Kenilworth Brisbane Lamington Glenugie lluka Kempsey 
Georissa laseroni 0 0 0 0 0 X 
Pleuropoma jana 0 0 0 0 X 
Pleuropoma draytonensis X X X 0 0 0 
Pupina wilcoxi X 0 0 0 0 0 
Pupina strangel X 0 0 0 0 0 
Pupinidae NN 1 0 0 0 X 0 0 
Pupinidae BR1 X 0 0 0 0 0 
Pupinella costata X 0 0 0 0 0 
Velepalaina strangel 0 0 X 0 0 0 
Tornatellinops jacksonensis 0 X 0 X X X 
Tornatellinops pressus 0 0 0 0 0 X 
Elasmias wekefieldiae 0 0 0 0 0 X 
Cylindrovertilla kingi 0 0 0 0 0 X 
Glyptopupoides egregia 0 0 0 0 0 X 
Gastrocopta pediculus 0 0 0 0 0 X 
Gastrocopta bifurcata X 0 0 0 X X 
Pupisoma circumlitum X 0 0 X 0 X 
Pupisoma porti X 0 X 0 0 X 
Pupillidae MV1 X X 0 0 0 X 
Pupillidae MV2 X 0 0 0 0 0 
*Lamellaxis clavulinus X 0 0 0 0 0 
*Lamellaxis gracilis X 0 0 0 X 0 
Strangesta ramsayi 0 X X 0 0 0 
Strangesta assimilans 0 0 0 X 0 0 
Strangesta bullacea 0 0 0 0 X 0 
Echotrida strangeoides X X X 0 0 0 
Saladelos macquariensis 0 0 0 0 0 X 
Saladelos urarensis 0 0 0 X 0 0 
Rhytididae MV1 0 0 0 0 0 X 
Rhytididae MV2 0 0 0 0 0 X 
Rhytididae MV3 X X X 0 0 X 
Rhytididae MV7 0 0 0 0 0 X 
Rhytididae NN1 0 0 0 X 0 0 
Rhytididae BR1 X X X 0 0 0 
Rhytididae SQ5 X X 0 0 0 0 
Hedleyella falconeri 0 X X 0 0 X 
Hedleyella maconelli X 0 0 0 0 0 
Brazieresta larreyi 0 0 0 0 0 X 
Pedinogyra rotabilis 0 X X 0 0 0 
lotula microcosmos X X X X X X 
*Paralaoma caputspinulae 0 0 0 0 0 X 
Punctidae MV1 0 0 X 0 0 0 
Ngairea corticicola 0 0 X 0 0 0 
Setomedea seticostata 0 X X 0 0 0 
Mussonula verax X X 0 0 0 0 
Nautiliropa omicron X X X 0 0 0 
Hedleyoconcha delta 0 0 X 0 0 X 
Gyrocochlea vinitincta 0 0 X 0 0 0 
Gyrocochlea convoluta 0 0 X 0 0 0 
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(continued) 


Kenilworth Brisbane Lamington Glenugie lluka Kempsey 


Gyrocochlea cinnamea 0 0 0 
Gyrocochlea planorbis 
Coenocharopa multiradiata 
Coenocharopa sordidus 
Coenocharopa parvicostata 
Coenocharopa macromphala 
Coenocharopa yessabahensis 
Rotacharopa densilamellata 
Allocharopa belli 
*Discocharopa aperta 
Egilomen globosa 
Egilomen lirata 

Elsothera genithecata 
Elsothera nautilodea 
Letomola contortus 
Rhophodon kempseyensis 
Rhophodon minutissimus 
Rhophodon colmani 
Omphaloropa varicosa 
species 1 indet. 
Charopidae BR2 
Charopidae BR11 
Charopidae BR13 
Charopidae BR15 
Charopidae BR22 
Charopidae BR23 
Charopidae BR28 
Charopidae BR29 
Charopidae BR30 
Charopidae BR31 
Charopidae BB39 
Charopidae BR41 
Charopidae BR42 
Charopidae BR43 
Charopidae NN6 
Charopidae NN7 
Charopidae NN8 
Charopidae NN17 
Charopidae MV 1 
Charopidae MV3 
Charopidae MV4 
Charopidae MV7 
Charopidae MV8 
Charopidae MV11 
Charopidae MV12 
Charopidae MV13 
Charopidae MV14 
Charopidae MV16 
Charopidae MV22 
Charopidae MV24 
Charopidae MV25 
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(continued) 


Charopidae MV26 
Charopidae MV29 
Charopidae MV30 
Charopidae MV31 
Charopidae MV35 
Charopidae MV46 
Charopidae NE3 
Charopidae SQ4 
Fastosarion aquila 
Parmavitrina planilabris 
Nitor medioximus 
Nitor pudibunda 

Nitor subrugata 

Nitor graftonensis 
Liardetia scandens 
Coneuplecta calculosa 
Tarocystis responsivus 
species 2 indet. 
Helicarionidae MV1 
Helicarionidae MV2 
Helicarionidae MV3 
Helicarionidae MV4 
Helicarionidae MV5 
Helicarionidae NN5 
Helicarionidae NN10 
Helicarionidae BR1 
Helicarionidae BR5 
Helicarionidae BR7 
Helicarionidae BR9 
*Zonitoides arboreus 
*Bradybaena similaris 
Sphaerospira fraseri 


Thersites novaehollandiae 


Thersites richmondiana 
Papuexul bidwilli 
Posorites turneri 
Posorites conscendens 
Ramogenia challengeri 


Austrochloritis separanda 


Austrochloritis porteri 


Austrochloritis nambucca 


Camaenidae SQ2 
Camaenidae BR1 
Camaenidae NN1 
Camaenidae NN2 
Camaenidae MV4 
Camaenidae MV6 
Camaenidae MV10 
Camaenidae MV11 


total 
unique 


% unique 


Kenilworth Brisbane Lamington Glenugie 


49 
20 
40.82 


STANISIC ET AL. 


0 


11.76 


0 
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43.9 


10 
58.82 


Iluka 


9 
45.45 


Kempsey 
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